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Abstract. Quarkonium production in proton-nucleus collisions is a powerful tool 
to disentangle cold nuclear matter effects. A model based on coherent energy loss 
is able to explain the available quarkonium suppression data in a broad range of 
rapidities, from fixed-target to collider energies, suggesting cold energy loss to be 
the dominant effect in quarkonium suppression in p-A collisions. This could be 
further tested in a high-energy fixed-target experiment using a proton or nucleus 
beam. The nuclear modification factors of J/r|) and T as a function of rapidity are 
computed in p-A collisions at -^5 = 114.6 GeV, and in p-Pb and Pb-Pb collisions 
at ^/s = 72 GeV. These center-of-mass energies correspond to the collision on 
fixed-target nuclei of 7 TeV protons and 2.76 TeV lead nuclei available at the 
LHC. 

1 Introduction 

Understanding the physical origin of quarkonium (J/i|), T) suppression in proton-nucleus (p-A) 
collisions has been a challenge for the past thirty years. This would of course be a prerequi¬ 
site in order to perform reliable baseline predictions in heavy-ion collisions, where quarkonia 
are expected to be dissociated due to Debye screening of the heavy-quark potential at finite 
temperature [1]. Perhaps even more importantly, the wealth of data (especially for J/t|) and T) 
available in p-A collisions could help to understand generic features of hard QCD processes in 
a nuclear environment. 

In everyday language, we often make the distinction between ‘fixed-target’ and ‘collider’ 
experiments when it comes to quarkonium production. This separation might look a bit artificial 
but not entirely: 

— In fixed-target experiments luminosities are often high, leading to abundant yields and 

consequently reduced statistical uncertainties. Moreover, thanks to the boost of the center- 
of-mass frame of the collision, the rapidity coverage of such experiments can extend up 
to very large values of rapidity (or Feynman-a;, Xp ~ x sinhy) using forward 

spectrometers. However, the highest energies ever reached are rather modest, = 38.7 GeV 
and ys = 41.6 GeV using respectively the 800 GeV and 920 GeV proton beams at the 
Tevatron and at HERA; 

— At collider energies - RHIC and LHG, to quote only the facilities accelerating heavy ions - 
unprecedented energies can be reached, respectively = 200 GeV and = 5.02 TeV, 
making for instance easier the production of T states, marginally measured in fixed-target 
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experiments. In terms of acceptance, quarkonia are detected in a narrow window in xp, 
centered around xp = 0. 


Let us illustrate this with a few examples, starting with one of the first experiments which 
measured J/i|) suppression in p-A collisions. The NAS spectrometer at the CERN SPS collected 
1.5 million J/ij) events (!) in hadron-nucleus collisions [2], allowing for precise measurements 
close to the kinematic edge of phase-space, xp ^ 0.75 (on the contrary, the coverage at RHIC 
and LHC is respectively [xpl < 0.2 and |xf| < 0.02 for J/rJj production). It is remarkable that 
these data, taken in the early 1980’s, prove as competitive as the most recent LHC results when 
it comes to understand J/ij; suppression in nuclei. More than a decade after NAS, the FNAL 
E866 experiment reported on high-statistics measurements of J/i|) (S x 10® events) and ij)' (10®) 
production on several nuclear targets in the range —0.2 < ccp < 0.9 [S]. These measurements^ are 
nicely supplemented by those carried out in d~Au collisions at RHIC (PHENIX [7,8], STAR [9]) 
and in p-Pb collisions at LHC (ALICE [10], LHCb [11]). At LHC, the relative suppression of 
T excited states (2S, SS) with respect to IS states has been performed by CMS [12], not to 
mention open heavy-flavor data (D mesons in ALICE [13], B mesons in CMS [14], non-prompt 
J/r[) coming from B decays in ALICE [15] and LHCb [11]). 

Several cold nuclear matter (CNM) effects could in principle alfect quarkonium yields in 
proton-nucleus collisions. Without being comprehensive, let us mention the following ones: 

— Quarkonia may interact inelastically with the surrounding nucleons they may encounter 
while propagating through the nucleus. Such nuclear absorption may happen when the 
quarkonium formation time (in the rest frame of the nucleus) is comparable or less than the 
medium length L, Tf x coshy^^^^ < L, where ry is the proper formation time (ry ~ 0.3 fm 
for both J/xj; and T), and is the quarkonium rapidity in the nucleus rest frame;^ 

— Parton distribution functions (PDF) are known to be different in a proton and in a nu¬ 
cleus at all values of x. Such nuclear PDF (nPDF) effects could either suppress or enhance 
quarkonium yields in p-A collisions (with respect to p-p collisions) depending on the value 
of the momentum fraction X 2 . When X 2 is small,® the nucleons in the nucleus act coherently 
leading to a reduction of the quarkonium yield in a nucleus - called shadowing [16], or 
saturation [17] to use a more modern language - as compared to the incoherent sum over 
A independent nucleons; 

— Nuclear transverse momentum broadening of the heavy quark pair induces coherent gluon 
radiation, arising from the interference between emission amplitudes off the initial projectile 
parton and the final color octet quark pair. This coherent medium-induced radiation leads 
to an average induced energy loss proportional to the quarkonium energy [18]. The conse¬ 
quences of coherent energy loss are quarkonium suppression (respectively, enhancement) at 
large positive (respectively, large negative) values of the rapidity and at all center-of-mass 
energies of the p-A collision. 

Obviously it is not because CNM effects can play a role that they necessarily do; in particular, 
the strength of each CNM effect is usually unknown a priori. A sound strategy is to investi¬ 
gate each of these effects separately, through a systematic and quantitative comparison to all 
available data, while keeping the smallest number of assumptions and free parameters. 

Quarkonium suppression reported at forward rapidities cannot be reproduced by either 
nuclear absorption or nPDF effects, nor by a mixture of both. Although the comparison to 
RHIC and LHC data only may still give the impression that strong nPDF effects could explain 
J/x]) data [10]^, the significant suppression measured by the fixed-target experiments (NA3 and 


^ as well as others measurements by, e.g., NASO [4], NA60 [5], and HERA-B [6], yet on a more 
restricted a;p range. 

® Note that is directly related to the momentum fraction X 2 carried by the target parton, 
coshx/i^jj = M±/{2mpX2), where M± is the quarkonium transverse mass and nip is the proton mass. 

® Typically when the time for the hard process to occur is large in the nucleus rest frame, th — 
(1/Mj_) X coshy,^j^ = l/(2mpa;2) ^ L. Using L = 10 fm, this would correspond to X 2 < 10“®. 

^ For examples of nPDF effects on quarkonium production in p-Pb collisions at LHC, see [?, 19]. 
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E866) on a wider xp range is clearly incompatible with the predictions of these two effects.® 
Without a doubt, the world data indicate that at least another cold nuclear matter effect is at 
play. 

Contrary to nuclear absorption or nPDF effects, the sole effect of coherent energy loss is able 
to reproduce the data on quarkonium suppression, from fixed-target to collider energies [22-24] . 
Detailed comparisons were published elsewhere, so let us only highlight the phenomenological 
successes of this approach: 


— The Xp (or, y) dependence of J/iJj suppression is well reproduced on a very large domain 
(up to large values of xp < 0.8, when data are available) and at all center-of-mass energies, 
from y/s ~ 20 GeV to y/s = 5 TeV; 

— The dependence is well reproduced too, either at a fixed-target experiment (E866) or 
at colliders (RHIC, LHC), although the dependence seems slightly more abrupt in the 
model than in collider data. The centrality dependence measured by PHENIX at RHIC is 
also nicely described; 

— T measurements in p-A collisions are compatible with the expected mass dependence of 
coherent energy loss, although the present experimental uncertainties are still fairly large; 

— Finally, an original prediction of coherent energy loss is a different magnitude of quarkonium 
suppression in p-A and tt-A collisions (in contrast with nuclear absorption effects, which 
should be independent of the projectile hadron), in agreement with the measurements of 
NA3. 


The strength of coherent energy loss depends on a single free parameter, namely the magnitude 
of the cold nuclear matter transport coefficient, Qq = 0.075 GeV^/fm at cc = 10“^, obtained 
from a fit of the precise E866 measurements in p-W collisions. 

We find it appealing that the variety of quarkonium measurements in p-A collisions can be 
described using a single CNM effect. Of course, by no means does this imply that no other CNM 
effects could play a role too, yet these clearly appear to be subleading when the quarkonium 
is produced at ‘large enough’ rapidity. Both nPDF and coherent energy loss effects could in 
principle be incorporated consistently in the picture. As a matter of fact, attempts have been 
made in [23]. However, given the large theoretical uncertainties on nuclear parton distributions® 
- especially for gluon densities at small x - we prefer to focus on the single (but in our opinion, 
leading) effect of coherent energy loss for which rather precise calculation can be performed. 

An exciting possibility to further constrain cold nuclear matter effects (on quarkonium 
production, but not only) would be to smash the LHC proton and lead beams on a collection 
of fixed nuclear targets [25]. We believe that this proposal would combine the above-discussed 
advantages of fixed-target and collider experiments. 

In this paper, the predictions for quarkonium suppression due to coherent energy loss in 
p-A collisions at ^/s = 114.6 GeV (corresponding to the nominal 7 TeV proton beam energy 
at the LHC) and Pb-A collisions at ^5 = 72 GeV (corresponding to the 2.76 TeV lead beam) 
are given. Before this, we recall in the next section the main ingredients of our approach. 


® An elegant way to be persuaded is to plot J/if) suppression data as a function of X 2 = Mx/y/s x 
exp(—j/) [20]. The suppression from either nuclear absorption or nPDF effects is expected to be a 
function of X 2 only, independent of ^/s, in violent disagreement with the accumulated data from fixed- 
target and RHIC experiments (see [3,21]). 

® This is due to the lack of small-x measurements in nuclear collisions. In this respect a high-energy 
electron-ion collider would be highly benehcial for the community. Let us mention in passing that 
no coherent energy loss effects are expected in deep-inelastic scattering experiments as the incoming 
projectile particle is color neutral [18]. 
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2 Coherent energy loss model in a nutshell 

2.1 Formulation 


We briefly detail in this section the basics of the model based on coherent energy loss used to 
describe (denoting J/r]) or T) suppression measured in proton-nucleus collisions.' The single 
differential p-A production cross section as a function of the ij) energy reads [23] 


AE 


{E) = r 

Jo 


, dcr^ 

deV{e,E,il)^{E + e) 


( 1 ) 


where E (respectively, e) is the energy (respectively, energy loss) of the QQ pair in the rest 
frame of the nucleus A. The upper limit on the energy loss is = mm{E,Ep — E), where 
Ep is the beam energy in that frame. V denotes the energy loss probability distribution, or 
quenching weight. 

The quenching weight is related to the medium-induced, coherent radiation spectrum dl/de 
given in [23] (and earlier in [18]), which is a very good approximation to the exact spectrum 
computed to all orders in the opacity expansion [27] . It depends on the accumulated transverse 
momentum transfer (assumed to satisfy <C M±) due to soft rescatterings in the 

nucleus, where L is the medium path-length and q the transport coefficient in cold nuclear 
matter. More precisely [23], 


q = qo 


10 


-2 


1 0.3 


min(a:o, X 2 ) 


xo = 


2mpL 


X2 = 




( 2 ) 


where y is the quarkonium rapidity in the center-of-mass frame of the proton-nucleon collision. 

Using the fact that the quenching weight is a scaling function of the variable e/E, namely 
E'P{e,E,t'^) = V{e/E,l'^), we can rewrite (1) as [26] 


1 do-pA 
A dy 


{y, v^) = y dSyV{e^^ - 1, q{y)L) {y + Sy, ^/s) . 


dcrl 


( 3 ) 


Here Sy'^^^{y) = min (ln2, j/^ax — y), with j/max = 'Al{^/s/M±) the maximal ij) rapidity (in the 
proton-nucleon c.m. frame) allowed by kinematics. Using (3) we can determine the nuclear 
modification factor in p-A collision. 


As mentioned in the introduction, quarkonium may suffer inelastic interaction with the 
surrounding nucleons while escaping the nucleus. Although we do not aim to include such an 
effect in the present predictions, we nevertheless indicate the critical rapidity 

(5) 

below which nuclear absorption might come into play. 

2.2 Ingredients 

The medium length L is obtained from a Glauber model calculation using realistic nuclear 
densities. The values are given in [23] and reproduced in Table 1 for the nuclei of interest in the 
present paper. In addition, Eq. (3) requires the knowledge of the p-p cross section. It is given by 

The model can also be formulated in heavy-ion (A-B) collision, see [26] for details. 
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a simple parameterization cosh yj , where the exponent n is obtained 

from a fit to p-p measurements. Lacking p-p data at the energies of interest = 114.6 GeV 
in p-A and y/s = 72 GeV in Pb-A collisions), an interpolation between the values obtained 
at FNAL {yds = 38.7 GeV) and RHIC (-y/s = 200 GeV) energies has been performed. The 
exponents used in the present paper are given in Table 2. Note that the normalization of the 
p-p cross section is irrelevant here as it cancels out when computing (4). 


Nucleus 

P 

Ga 

Gu 

Pb 

Atomic mass 

1 

40 

63 

208 

L (fm) 

1.5 

5.69 

6.67 

10.11 


Table 1. Values of L used in p, Ca, Cu, and Pb targets. 


Mode 


Pb-A p-A 


^ (GeV) 72 114.6 

n,,, 5.1 ±0.2 6.0 ±0.3 

J/4> 

4.1 ±0.3 5.0 ±0.4 


Table 2. Values of n used at ^/s = 72 GeV and ^/s = 114.6 GeV for 1/4) and T. 


The transport coefficient % is the only free parameter of the model. It is determined by 
fitting the J/r)) suppression measured by E866 [3] in p-W over p-Be collisions {^/s = 38.7 GeV), 
see [23]. The obtained value is go = 0-075lg;Qj5 GeV^/fm. 


3 Results 

3.1 p-A mode 

The predictions for J/i|) and T suppression in p-Ca, p-Cu and p-Pb collisions at y/s = 
114.6 GeV are shown in Figure 1. The rapidity range is chosen to match the acceptance of 
detectors like LHGb. In terms of Feynman-x, the range —3 < j/ < 1 (respectively, —2<y<l) 
correspond to —0.54 < a;F < 0.06 (respectively, —0.60 < ccp < 0.19) for (respectively, T). 

The J/r[> suppression is rather moderate, less than 20%, and does not vary too strongly 
with rapidity except at very negative rapidity values, i/ < yg — ~2, where J/ij) enhancement 
(i?pA > 1) can be seen. In this rapidity region, however, nuclear absorption may come into 
play as can be seen from the vertical arrows indicating the values of — —1) for each 

target. 

The shape of T suppression is similar. The value of the rapidity at which i?pA(2/o) = 1 is yo — 
— 1, j.e., one more unit than in the J/i|) channel. This can be understood from the approximate 
xf scaling present in the model [23],® which would predict the difference between these two 
‘crossing points’ to be yQ — yg^’^ ln(M^/M'^/^) ~ 1.1. Once more, nuclear absorption may 

affect T suppression, although maybe not as much as the J/r|^ because of its smaller radius. 

® At a given rapidity y, the corresponding value of xp is larger for T than for J/4> due to the larger 
transverse mass, a^p oc M±. 
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Fig. 1. J/'i]) (left) and T (right) suppression in p-Ca, p-Cu and p-Pb collisions at = 114.6 GeV. 


3.2 A-p mode 

Let us move now to calculations corresponding to an incoming 2.76 TeV Pb beam on a pro¬ 
ton and a Pb target, shown in Figure 2. This configuration allows for probing more easily 
quarkonium suppression in the proton fragmentation, i.e., at positive rapidities.® The chosen 
rapidity range —2<y<2 (respectively, —1.5 < y < 1.5) correspond to —0.31 < a:F < 0.31 
(—0.56 < Xp < 0.56) for J/t)j (T). The lower center-of-mass energy however shifts y™* in p-Pb 
collisions towards larger values, possibly leading to more pronounced nuclear absorption. 

In Pb-Pb collisions the suppression is naturally an even function of y. In such collisions, 
one expects a hot medium to be formed leading to extra quarkonium suppression. Therefore 
the results should rather be seen as baseline calculations than genuine predictions. Moreover, 
in A-A collisions the condition for hadronization taking place outside both nuclei reads < 
y < — . This condition is only met in Pb-Pb collisions at ^/s = 72 GeV around mid-rapidity, 

\y\ < 0.1. At larger \y\, the quarkonium state shall be fully formed in one of the two nuclei, and 
thus possibly sensitive to nuclear absorption. 




Fig. 2. J/\l) (left) and T (right) suppression in p-Pb and Pb-Pb collisions at = 72 GeV. 


® Although the Pb nucleus collides on a proton, we shall keep the convention that positive values of 
y correspond to the proton fragmentation region. 
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